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We calculate nucleon-nucleon scattering at low energies and large impact parameter 



in the Skyrme model within the framework for soliton scattering proposed by Manton. 

ON 

This corresponds to a truncation of the degrees of freedom to the twelve collective 

<^ . coordinates which essentially describe the rigid body motion of the pair of Skyrmions. 

We take to its logical conclusion the result that the induced kinetic energy for these 

Ph, collective coordinates in the product ansatz behaves as one over the separation and 

Ph| 

D , hence can dominate over the potential. This behaviour implies to leading order that 

we can drop the potential and the resulting motion reduces simply to geodesic motion 
on the manifold parametrized by the variables of the product ansatz. We formulate 
the semi-classical quantization of these variables to obtain the motion corresponding 
to the nucleonic states of the Skyrme model. This is the appropriate description 
for the nucleons in order to consider their scattering within Manton's framework in 
the semi-classical approximation. We investigate the implications for the scattering 
of nucleons with various initial polarizations using the approximation method of 
"variation of constants". 
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1. Introduction 

The Skyrme mo de£l is unique among models describing the strong interactions. It is con- 
ceived from the very basic and fundamental principles concerning the symmetries and their 
realizations in the strong interaction, and it can be used to make quantitative predictions 
about low energy hadronic dynamics. It contains exactly two parameters, f n and e, which 
can in principle be computed from QCD. In practice, they are fitted by phenomenological 
consideration. Once these parameters are fixed, the model can make reasonably accurate 
predictions in the meson sector§ and for the static properties of single baryonsi. 

The case B = 2 is very interesting since it contains both the scattering of two nucleons 
and their bound state, the deuteron. Much effort has been invested in extracting the nucleon- 
nucleon potentiali. This comprises of a calculation of the potential energy of two well 
separated Skyrmions, taken in the product ansatz or slight variations thereof, completed with 
a projection onto asymptotic quantum nucleonic states. For large separation d between the 
nucleons, the static potential behaves as 1/d 3 , for the case of massless pions, and contributes 
only to the spin-spin and tensor channels of the nucleon-nucleon potential. Adding a pion 
mass cuts off all interactions with the usual Yukawa decaying exponential, however the 
potential now has a leading 1/d behaviour in the tensor channel and a 1/d 2 central channel, 
albeit, both are multiplied by two powers of the pion mass. The product ansatz is deficient 
to produce an intermediate range attraction in the central channel, however, it has been 
demonstrated by exact numerical calculations unfettered by this ansatz, that the Skyrme 
model does indeed contain a central attraction!'!. The Skyrme model has not been too useful 
for determining the potential at small separation. The exact deuteron profile is thought to be 
a toroidal configuration where the Skyrmions loose their identity. In this way it is clear that 
the product ansatz will fail to properly describe the situation, since it distinctly identifies 
the two Skyrmions. 

It is the aim of this paper to consider the influence of the kinetic energy on the nucleon- 
nucleon interaction. Within the framework of MantonB for soliton scattering, the system is 
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truncated to the finite number of degrees of freedom parametrizing the low energy critical 
points of the static configuration and the gradient flow curves joining them together. The 
full potential serves to induce a potential on this sub-manifold, while the full kinetic energy 
serves to induce a metric. We, and also Schroer s0 have computed the induced kinetic energy 
for the case of massless pionsi and found a leading term which behaves as 1/d. This term 
can dominate for large separation over the contribution from the potential which behaves as 
1/d 3 . We wish to investigate the effect of this kinetic term on nucleon-nucleon scattering. 
There are two questions to address: firstly how to properly formulate the inclusion of the 
term within the low energy, semi-classical framework and secondly, how to characterize 
its effects on the scattering. We will work with the hypothesis that the potential term 
is in fact negligible, we discuss in the conclusion the actual physical realization of this 
eventuality. With this hypothesis, the scattering corresponds to geodesic motion in the 
induced metric on the sub-manifold. The effect of this induced metric is to modify the 
equations of motion in a rather complicated way. The utility of the differential geometric 
formulation is less evident with respect to the corresponding work on scattering of BPS 
monopolesi or self dual vortices^, primarily because of the lack of symmetries and useful 
algebraic structures. Instead the equations are better formulated as one complicated, exact, 
dynamical equation for the coordinate corresponding to the relative separation, augmented 
by conservation laws, which are now violated by perturbative corrections coming from the 
new interaction, for the generators of rotation and for isorotation of the free system. We 
make a further approximation, using the method of "variation of constants" to compute the 
scattering. The method essentially consists of using the free trajectories in the interaction 
terms to produce evolution of the previously conserved (constant) quantities. We find non- 
trivial scattering for all initial polarizations. 

The extraction of the implied nucleon-nucleon scattering is non-standard. Up to now, 
most investigations have focussed directly on the projection onto asymptotic, quantum nu- 
cleonic states. This requires a consistent interpretation of the kinetic energy as supplying 
an induced metric on the sub-manifold. This metric modifies the generally covariant Lapla- 
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cian denned on the sub-manifold which then should be included in the relevant Schrodinger 
operator!!. We find that this procedure has not been fully implementedi'0'El. It remains 
to be seen whether directly quantizing the truncated system is phenomenologically viable. 
Actually we find this procedure somewhat orthogonal to the incorporation of classical soli- 
tons within the quantum theory. This incorporation is necessarily formulated in terms of the 
semi-classical perturbative expansion, indeed the soliton is a classical solution. In this article 
we give the construction of the semi-classical nucleonic states using the Bohr-Sommerfeld 
(WKB) type quantization rules. These states are then amenable to calculation of the scat- 
tering within the semi-classical approximation. 

In section 2 we give the equations of motion and the approximation methods used. In 
section 3 we give the analysis of the semi-classical quantum nucleon states. In section 4 we 
summarize our results for the scattering and we terminate with our conclusions in section 5. 



2. Equations of Motion and Approximation Methods 

In this section we will indicate briefly the path to deriving the equations of motion and the 
approximations that we use in our treatement. For a more detailed description we refer the 
reader to our previous article and references thereini. 
With the definitions 

CUV) = ~\ trnurf] (1) 

n;(u) = -- tr[r a u^d,u} (2) 

D ab (U) = \ tr[T a Ur b tf] (3) 

where the field U (x, t) is an element of the 2x2 matrix valued representation of SU(2), the 
Skyrme Lagrangean becomes 

C sk = f -fC,{U) ■ C(U) - \c,{U) ■ C»{U)C V {U) ■ C U {U) - Cp(U) ■ C{U)C V {U) ■ C(U) . 

(4) 



We replace for U(x, t) in this Lagrangean with the product ansatz 

U(x, t) = AU s (x - RjA^BUsix - R 2 )B^ (5) 

where U s (x) = e if ^- f is the single Skyrmion with profile function f(r), A(t),B(t) G SU(2) 
define the respective spin and isospin orientations of the Skyrmions and R\ and R 2 are their 
positions. This corresponds to the truncation^ to the relevant degrees of freedom for the low 
energy behaviour. We obtain the kinetic energy by keeping only the terms which contain 
time derivatives, these are automatically quadratic, the remaining terms contain no time 
derivatives and comprise the potential. 

The product ansatz makes sense only for well separated Skyrmions, when they are close, 
we expect strong deformations and, in fact, the minimal critical point is conceived to be a 
toroidal configuration^. A proposal to globally, approximately describe the sub-manifold of 
low energy critical points and the union of gradient flow curves joining these together has 
been put forward by Atiyah and MantonS We do not consider, however, the possibility 
of close encounters between the Skyrmions and hence the product ansatz should yield at 
least the correct, leading asymptotic result. The induced kinetic energy has a well defined 
expansion in inverse powers of the separation d. The zero order terms correspond to the 
free kinetic energy of each individual Skyrmion. The leading correction behaves as 1/d 
and contributes to a spin-spin and tensor interaction between the Skyrmions, as we! and 
SchroersB have found. Thus the kinetic energy is given by 

T = \Md 2 + 2k(C a (A) C a (A) + C a {B) £ a (B)) 

( 6 ) 

+ A e iac e jbd ftd( B ) ( S a _ d i d jj D ab (A^B) + 0(1/ d 2 ) 

where M is the Skyrmion mass, A its moment of inertia, A = 2iiK 2 f 2 , F(r) ~ n/r 2 at large 
r, C a (A) = C*{A), TZ a {A) = 111(A) and d = d/d = (R x - R 2 )/d. The potential has been 
calculated by many authorsi'lllii. It has a leading term which contributes to the tensor 
channel and behaves as 1/d 3 , apart from the "free" contribution which gives the mass of 
each Skyrmion: 
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V = 2M + 4vr/> 2 



(l-cosfl)(3(n-d) 2 -l) 



(7) 



The complete Lagrangean of the system is given by 



L 



T - V. 



(8) 



In this paper we focus on the effects of the kinetic term. The equations of motion arising 
from the kinetic term are neither illuminating nor practical. It is useful just to write down 
the approximate equations. The free Lagrangean is completely integrable and corresponds 
exactly to two free, spherically symmetric tops. We can express the approximate equations 
as corrections to conservation laws using the approximation of Lagrang Here we take 
the Poisson brakets implied by the free dynamics to compute the corrections implied by the 
interaction to the equations of motion. This yields a practical set of equations which must 
be treated, however, with further perturbation measures 



f t C k {A) = ^e iac ei bd K c {A)K d {B){5 i: > - d i d j )e ke fD fa (A)D eb {B) 

iC k (B) = ^e^ 6d ft c (A)ft d (5)(<^ - <M)e fcei 'D ae (A^)D fb (B) 

f t n k (A) = -^e iac e jbd K d (B)(8 i i - tft£) e kcf K j \A)D ab (A^ B) + e ka W fb (A^B)n c (A) 

f TZ k (B) = -^ 2 e iac ^ bd TZ c (A)(5 i ^ - $d?) \e kd fn^{B)D ah {A^B) + e kb f D af (A^B)TZ d (B) . 



The method of variation of constants is ideal for the following analysis. Here we replace 
the variables in the right hand side of the "conservation equations" by their trajectories 
coming from the free dynamics. This gives rise to a "variation" of the previously conserved 
"constants". The procedure can be iterated to give higher order corrections. One should 
however always maintain consistency with the first (Lagrange) approximation. In fact, the 
method of "variation of constants" is only useful with respect to the equation for d, where it 
gives the scattering trajectory. The change in the spin or isospin governed by the subsequent 
four equations cannot be used in conjunction with this approximation method because of the 
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5 ij d k + 6i k <P + 5 ik d j - 3dWd k e i<w e jbd TZ c (A)TZ d (B)D ab (A^B) 



(9) 
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long range of the interaction. The results give an infinite change in these angular momenta, 
which is not reliable. 



3. Semiclassical Nucleon states 

In this section we will construct quantum states corresponding to the nucleons, using the 
semi-classical Bohr-Sommerfeld quantization rules applied to the classical motion of the free 
Skyrmion. To apply these rules, we must find the classical action-angle variables of the 
system and impose quantization conditionsE!!!. This enables us to treat a Skyrmion as 
a nucleon consistently within the context of the semi-classical expansion. Going directly 
to the quantization of the SkyrmionlHIS is perhaps not in the spirit of the semi-classical 
formulation which is, we believe, incumbent upon us once we have begun to consider classical 
solutions. In addition, we find that with the proper semi-classical quantization we can make 
physical sense of the specific, classical Skyrmion scattering trajectories for the scattering of 
nucleons. We start with the Lagrangean for a single Skyrmion in the center of mass frame 
which is given by: 



where d« is the time derivative of the eij from the parametrization A = a + id ■ r with 
Oq + a 2 = 1 of the matrix defining the rotational characteristics of the Skyrmion. Following 
the treatment of the motion of a rigid body in space, we go to the system of coordinates 
corresponding to the usual Euler angles: 



where a, 7 G [0, 2n] and (3 G [0, 7r]. The a, are related to the Euler angles by the equations: 



L = -M + Air[AU] 
= -M + 2\ E 3 i=0 d 2 



(10) 



A = e - iar 3/ 2 g-«/3T2/2 e _i 7 T 3 /2 



(11) 
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a = cos(§) cos(^) 
ai = sin(f ) sini 9 ^) 
a 2 = — sin(|) cos( 9 ^f L ) 
a 3 = -cos(f)sin(2p). 

Replacing them into the Lagrangean ([[!]) we get 

L = — M + ^A [a 2 + /? 2 + 7 2 + 2«7 cos (3} (13) 

which is the well known Lagrangean for rigid body motion. In this case, the Euler angles 
are used to describe the motion of the body fixed axes relative to the laboratory fixed axes. 
In our case, the identification of the Euler angles will be completely different. 
The action-angle variables Jj are defined by0 



J; 



jpidqi (14) 



where Pi is the momentum conjugate to the coordinate g«, and the integral is taken along a 
closed path followed by the system during one period in the plane (qi,Pi) of phase space. In 
our case there are three such variables, one for each angle. Due to the cyclic nature of the 
angles a and 7, J a and J 7 are readily computed: 

J a — j> ' Padq a = 2np a = 27rA[a; + 7COS/3] (15) 

J 7 = j> p^dq^ = 2npj = 27rA[7 + dcos/3]. (16) 

In contrast, Jp is much more complicated for general orbits, however, for the type of motions 
that we will consider it is zero. The motion of a free Skyrmion corresponds simply to uniform 
isorotation about a fixed axis. The equations of motion of the Lagrangean flT3| ) are 

a + 7 cos (3 — 7/3 sin (3 = 

(3 + 2d 7 = (17) 
7 + a cos (3 — a (3 sin (3 = 



Without loss of generality we choose the axis of isorotation to be the third axis in isospace. 
Then the time evolution is clearly given by linear time dependence of a and/or 7 with /3 — 

Or 7T. 

The Bohr-Sommerfeld quantization condition, derived by path integral methods@@ is 
W = J2 J i = ( n + OK n = 0,±1,±2,---. (18) 

i 

£ is a correction factor arising from the functional integral over Gaussian fluctuations about 
the classical trajectory. In zero ordei@ we can neglect £. Then we find 

J a + J^ = nh, n = 0,±1,±2,--- (19) 

in perfect accord with the original Bohr-Sommerfeld quantization rules. In order to precisely 
fix the time dependence of a and 7 corresponding to nucleon states we must make reference 
to the isospin and spin generators. These are given by, respectively: 

£3 = 2A£ 3 = 

277 (20) 

n 3 = 2\n 3 = 

Z7T 

and satisfy the usual angular momentum algebra. The symmetries of the Skyrmion dictate 
that spin is equal to isospin, thus 

\Ja\ = \J-y\- (21) 

Then applying the quantization condition with n = 1, gives 

|£ 3 | = |£ 3 | = - (22) 
We will follow the convention of Adkins et al i who define the isospin generator as 
h = -£ 3 (23) 
and the spin generator as 

J 3 = n 3 . (24) 



Consider first the case /3 = 0. Replacing this in equation (|15D , (|16|), (pOl), ( p3|) and (^4j) we 

get 

I 3 = A(d + 7) = 2Xu = ±1/2 

V ^ ' (25) 

J 3 = -A(d + 7) = -2\u = =pl/2 

with (a + 7)/2 = <fi(t) = ut + O - This type of angular motion produces Skyrmions with 
spin and isospin antiparallel, and eventually the states \p i> and \n |>. The proton state 
corresponds to isospin +1/2 along the 3 axis in isospace while the neutron corresponds to 
isospin —1/2. The selection of the third axis, we emphasize, is purely convention. States 
polarized along some other axis in isospace simply corresponds to linear superpositions of 
our states. The corresponding beams would be a mixture of protons and neutrons. 
Replacing this solution in equation ( |12"1) gives 



A = cos (bit) — i sin d(t)T 3 

W W (26) 

_ e -i0(t)r 3 /2 

where uj > corresponds to the state \p |> and u < to \n f >. 
For (3 = 7T 

7 3 = A(d - 7) = 2Ac^ = ±1/2 

1 V ' (27) 

J 3 = -A( 7 - a) = 2Xdj = ±1/2 

with (a — 7)/2 = if)(t) = ut + ipo and the corresponding matrix 
A = ilsmip^r 1 — cos-^(t)T 2 ] 



(28) 



-ii/j(t)T 3 /2 2 „#(t)r 3 /2 



represents the state \p |> with uj > and \n [> with uj < 0. We can see the similarity of the 
forms for the matrices A corresponding to nucleon states and the wave functions obtained in 
the directly quantum version of Adkins et al i. In the next section we study the scattering 
of the nucleons as described by these states, using the geodesies approximation. 
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4. Nucleon-Nucleon scattering 

The nucleoli states found in the previous section are ideally suited for our treatment of 
nucleon-nucleon scattering using the approximation method of the "variation of constants" . 
We simply replace into the right hand side of equation (^) the quantized classical trajectories 
found in the previous section. To calculate the change induced in the previously constant 
quantities we integrate the equations from t = — oo to t — +00 over one free trajectory This 
is useful for the equation governing the relative motion (|^). For the angular variables the 
perturbation method is too crude. It is perhaps a reasonable perturbative scheme locally in 
time, however, to integrate the motion over all time leads to infinite changes in the various 
angular momenta. We do not trust such results and we suggest that the approximation 
should be iterated several times to obtain a better approximation to the local (in time) 
perturbation of the free trajectories. This could then perhaps be integrated over all time. 
At the root of the divergence is the long range (1/cf) nature of the interaction term. Clearly 
for a physical theory with non-zero pion mass, this interaction would be cut-off by the usual 
Yukawa factor, rendering all such integrations finite. Thus we will consider the variation of 
the angular momenta only in so far as making the observation that the equations of motion 
with the interaction imply non-trivial spin flip and charge flip scattering. These effects are 
compatible with the exchange of charge carrying (pions) and spin carrying (vector mesons) 
intermediate particles. The first equation of @ which describes the relative position is of 
course sufficient to obtain the exclusive scattering. We present below the results for the 
scattering of specially polarized nucleons, described in our semi-classical formalism. 

It is interesting to observe that in our formalism an additional parameter arises which 
describes the initial state of two incoming, polarized nucleons. This parameter, along with 
the impact parameter, the initial velocity and the direction of polarization, actually selects 
the particular scattering trajectory followed by the nucleons. The parameter describes the 
relative orientation of the Skyrmions at a fixed time. It plays in fact a role similar to a 
hidden variable. An incoming pair of physical nucleons, in our formalism, has a fixed value 
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for this parameter, which is only "measured" after the scattering takes place. In a physical 
experiment consisting of incoming beams of nucleons giving rise to collisions or scattering 
of pairs of nucleons, the value of this parameter will be uniformly distributed. A similar 
parameter arises in the case of the scattering of BPS monopoles0. 

There is an immediate separation of the scattering into two cases, depending on whether 
D a b(A^B) is time independent or not. When it depends on time for large values of the ratio 
(uj/v) there is an exponential suppression of the scattering where v is the relative velocity 
and 7 the impact parameter. This is quite evident, for slowly translating Skyrmions, the 
prescribed rotations imposed by selecting semi-classically quantized nucleon states have the 
effect of averaging the interaction to zero. 

We find the following expression for time variations of the previously constant relative 
momentum p = (M/2)d: 



D a b(A !i B) time independent: 



V \P T 



n l n I d k Au 2 1k 

-P k = — cos{25)d k 29 

I I dt GT 

P IP I 

n t n t 



Ap k = ~^{d k + Ar k f ■ d - 6d k (f ■ df] 
P T V I " " (30) 



sin(<5), cos(5), 0) 



Apk = -^[d k + Ar k r ■ d - 6d k (r ■ d) 2 } 
n t n I ' (31) 

r k = (sin(<5), - cos(<5), 0) 



D a b(A^B) time dependent: 

p ^ n ^ d k Au 2 , A ^ /00 , 
— p = cosUujt + e)d k (32) 

i * dt d 1 
pin] 
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f t pk = ^\d k + 4r k f -d-6d k {f- d) 2 } 
P T n T ' (33) 

r k = (- sin(2cjt + e), cos(2wt + e), 0) 



f p k = ^ d k + Ar k f .d-6d k (f- df] 
V in i ' (34) 

r k = (- sin(2cut + e), - cos(2wt + e), 0) 

with d(t) in the right hand side given by 

(d x ,d y ,d z ) = (ui,7,0). (35) 

A(i) and -B(t) are given by the angular dependence found in the previous section corre- 
sponding to the associated nucleon states, with 5 = 0q — 0^ and e = 0q + 0^ in self evident 
notation. The right hand sides can be interpreted, at this level of our approximation, as 
coming from a spin-spin channel and a tensor channel interaction. We stress that this is only 
a correspondance, the true effect of the kinetic term is to supply a non-trivial connection in 
the geodesic equations on the low energy sub-manifold and not to modify the potential. 

To find the actual change mp and hence the scattering angle, we integrate these equations 
from t = — oo to t — +oo. The expressions for the scattering of p on n each contain a time 
dependent A^B. When integrated these yield an exponentially suppressed variation in the 
dimensionless group (wy/v) 

~ e~^\ (36) 

Thus in the limit v — > we get negligible scattering in these cases. For the cases pp or 
nn we find the scattering angle depends on a variable 5 which corresponds to the phase lag 
between the rotation of A(t) and B(t). 

We give the p | p | scattering in detail to make things concrete. Here we find 

p\t) = -A cos(25) { jZ%> dt + f vk (3?) 

= -A UJ 2 C0S(25) (-„ 2( „2t2" +7 2)l/2 + 7 2(„2 t 2+ 7 2)l/2) + 

This yields 
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p k (+oo) = -2A to 2 cos 25 A- + — v k (38) 

7 2 t> 2 



from which we calculate the cosine of the scattering angle 

(39) 



cos 9 = p(+ 00 )-p(~ °) 

|p(+oo)||p(-oo) 



M~/v 2 1 

This is the first analytical calculation of nucleon-nucleon scattering from essentially first 
principles, without recourse to ad hoc models or potentials. We reemphasize that the Skyrme 
model is in principle derivable from QCD and /„- and e are, as such, calculable parameters 
and, in that sense this is also a QCD calculation. 

We wish to point out that in the limit that the initial velocity vanishes, for fixed uj 
and 7, we recover 90° scattering. This is, however, not so surprising as it is a property 
also shared by the Coulomb and many other interactions treated with our approximation. 
90° scattering is hardly remarkable except at zero impact parameter, where of course, it is 
impossible to avoid the region of close proximity of the nucleons and it seems important 
that the configurations pass through the minimal, toroidal configuration. 

We have up to now considered scattering in the (x, y) plane with spin polarized in the 
orthogonal z direction. The tensorial nature of the interacton implies that the forces depend 
on the angle between the axis of separation and spin polarization. If we choose the spin 
polarization along an axis tilted with respect to the normal to the initial scattering plane 
we get complicated, three dimensional scattering trajectories. 



5. Conclusions 

In this paper we have taken to its logical conclusion the study of nucleon-nucleon scattering 
in the Skyrme model using the method of truncation to the relevant finite number of degrees 
of freedom and their ensuing dynamics as implied by an induced metric and potential. We 
have restricted ourselves to large separations between the Skyrmions. In this regime, the 
Skyrmions are well described to leading order by the product ansatz and the manifold of 
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collective coordinates is parametrized by the variables of the product ansatz. The induced 
metric, calculated to leading order behaves as 1/d while the induced potential behaves as 
1/d?. In principle there can be a region where the induced metric contribution dominates 
and we can neglect the potential. It is exactly this point of view which we have taken to its 
logical conclusion. We find that the metric induces an interaction which can be interpreted, 
within our approximation method, as a spin-spin and a tensor interaction. Unfortunately 
it seems that the domination by the metric term is not physically realized. The induced 
kinetic term is multiplied by essentially the frequencies of angular rotation of the Skyrmions 
while the potential term has two extra powers of the separation in the denominator. Thus 
for the kinetic term to dominate, the frequencies should be much larger than the separation, 
in units where distance, time and energies are all measured in MeV. The frequencies are 
fixed by the Bohr-Sommerfeld quantization rules which give 

u ~ lOOMeV. (40) 

This corresponds to a region of validity for a separation of about 3 fm and greater. However 
there is much latitude available since the values of f n and e which go into determining oj are 
fixed only by choosing two experimental inputs. f n , e can vary as much as 10-30% thus we do 
not feel overly concerned with exceeding the regime of validity. Our approximation would of 
course be better justified for the case of A— A scattering where u ~ 300MeV corresponding 
to a separation of 1 fm. In any case, we do believe however, that it is not physically 
reasonable to consider the scattering of nucleons with the metric term alone and we expect 
a contribution from the potential term which is of the same order of magnitude. We do not 
expect, however, any great, qualitative modification of the scattering upon inclusion of the 
potential term, it is of similar strength but actually contributes only in the tensor channel 
for the case of massless pions. 

We find our work as a prescriptive and requisite endeavour to investigate the nucleon- 
nucleon interaction arising from the kinetic term. We have shown how to properly formulate 
the nucleon states within the semi-classical approximation. We have treated the scattering 
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and computed the scattering angles in a systematic perturbative method. Future work 
should include consideration of a non-zero pion mass, which leads to a central channel 
interaction, a better control of the perturbation method, a departure from the product ansatz 
and a proper treatment of the region of close proximity, to test the validity of our formalism 
in the phenomenology of low energy nucleon-nucleon scattering and of the static quantum 
states in this sector. 
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